Magnetic vortex lattice in HgBa2Cu04+5 observed by small-angle neutron scattering 
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We report a direct observation of the magnetic vortex lattice in the model high-temperature 
superconductor HgBa2Cu04+i. Using small-angle neutron scattering on high-quality crystals, we 
observe two equal domains of undistorted triangular vortex lattices well-aligned with the tetragonal 
crystallographic axes. The signal decreases rapidly with increasing magnetic field and vanishes above 
0.4Tesla, which we attribute to a crossover from a three-dimensional to a two-dimensional vortex 
system, similar to previous results for the more anisotropic compound Bi2.i5Sri,95CaCu208+a • Our 
result indicates that a triangular vortex lattice at low magnetic fields is a generic property to cuprates 
with critical temperatures above 80 K. 

PACS numbers: 74.25.Uv, 61.05.fg, 74.72.Gh 



I. INTRODUCTION 

The cuprate high-temperature superconductors ex- 
hibit a rich variety of mesoscopic phenomena associated 
with magnetic vortices in the mixed state>i In clean sam- 
ples, the vortices lower their energy by forming a long- 
range-ordered vortex lattice (VL), which provides an op- 
portunity to study the superconductivity of the bulk. In 
a simplified picturefS, the vortex physics is determined 
by two parameters: the penetration depth A of mag- 
netic field into the superconducting phase and the co- 
herence length ^ of the superconducting order parame- 
ter. Since both of these parameters exhibit pronounced 
anisotropies in the cuprates, the VL is expected to reflect 
important characteristics of the superconductivity, such 
as the d-wave pairing symmetry, ^^'^ However, there exists 
a substantial variety of VL properties among different 
compounds, including Bi2.i5Sri.95CaCu2p8-i-5 (Bi2212, 
Ref. 0), YBaaCugOg+s (YBCO, Refs. i and i), and 
La2-^Sr^Cu04 (LSCO, Ref. [13). The situation is fur- 
ther complicated by the presence of twin domains and 
boundaries in the most studied compound YBCO.— 
Therefore, it is highly desirable to measure a structurally- 
simple representative system in order to distinguish the 
generic vortex physics from material-specific aspects. 

We report here the first observation of the VL 
in HgBa2Cu04+5 (Hgl201) using small-angle neutron 
scattering (SANS). IIgl201 is in many ways a model 
systemjiiii^ First, Hgl201 exhibits the highest supercon- 
ducting critical temperature (Tc) among all cuprates with 
a single Cu02 layer per primitive unit cell. Notably, our 
samples exhibit the highest Tc among all materials in 
which VL studies by SANS have been carried out to date. 
Second, the Cu02 layer of Hgl201 is free of long-range 
structural distortionSfi^ and disorder resides in the Hg- 
O layer, relatively far away from the Cu02 layer.—"— 
As disorder and defects are common causes for vortex 
pinning, their relative absence in IIgl201 results in weak 



pinning - and is important for the formation of the VL. 
Third, Hgl201 bridges some differences between other 
cuprates: LSCO possesses a single-layer structure and 
relatively low Tc, whereas YBCO and Bi2212 are double- 
layer compounds with a high Tc. On the other hand, 
Hgl201 possesses a high Tc and it is a single-layer com- 
pound. Finally, while the other three compounds exhibit 
long-range structural distortions in the superconducting 
doping range, the high tetragonal symmetry of Hgl201 
(and the absence of long-range distortion) allows us to 
avoid complications arising from vortex-pinning due to 
crystallographic twin boundaries. 



II. EXPERIMENTAL METHODS 

High-quality Hgl201 crystals have recently become 
available and proven to be suitable for cutting-edge neu- 
tron scattering measurements f22r— The current study 
was carried out on two arrays (total mass of approxi- 
mately 60 and 90 nig) of nearly optimally-doped Hgl201 
crystals. The crystals were grown by a self-flux method?^ 
and subsequently heat-treatedi^ for up to 180 days in or- 
der to reach the desired homogeneous hole concentration. 
They were further selected according to the results of 
magnetometery measurements to ensure high qualityji^ 
and had a uniform Tc of 94 K (determined from the tran- 
sition mid-point) with a width of less than 2 K. The crys- 
tals were co-aligned and glued on a silicon wafer using 
an x-ray Laue diffractometer. The c-axes of the crys- 
tals were within 1° perpendicular to the surface of the 
wafer, and the a- and 6-axes were co-aligned to within 
1.5°. Since the two samples gave essentially the same re- 
sult, we do not distinguish between them in the remain- 
der of this paper. The measurements were performed 
on the SANS-I spectrometer at the Paul Scherrer Insti- 
tute, Switzerland. Incident neutrons of wavelengths (An) 
of 5-10 A were selected and coUimated over a distance 
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of 8-18 m before the sample. Diffracted neutrons were 
observed by a two-dimensional multidetector located 13- 
20 m after the sample. The sample was mounted in a 
horizontal cryomagnet with the a-axis vertical and the 
c-axis pointing along the direction of the applied mag- 
netic field (within 1°). The VL was prepared by cooling 
the sample from above Tc in a field oscillating around 
the target value (by ~ 1% of the target value). This 
method^ improved the VL ordering by keeping the vor- 
tices away from local pinning potentials. After reaching 
the lowest temperature of 2K, the field was held sta- 
tionary during diffraction measurements, which involved 
rotating the sample and cryomagnet together to angles 
that brought various diffraction peaks onto the Bragg 
condition at the detector. Measurements at tempera- 
tures between 2 K and Tc were performed upon heating 
without changing the field. For all the data presented 
here, the background measured at T = 2 K in zero field 
or at T > Tc has been subtracted from the field-cooled 
measurements to extract the VL signal. 



III. RESULTS 

Figure [T] shows diffraction patterns collected at T = 
2K for B — 0.2 and 0.3 T. The data are averaged over 
multiple sample rotation angles to display all first-order 
diffraction peaks. Given the tetragonal crystallographic 
symmetry of Hgl201, the VL is expected to form two 
equally-populated domains aligned with the a- and b- 
axes. The presence of twelve (2 x 6) diffraction peaks with 
comparable intensities confirms this, and indicates that 
the VL is triangular. Well-defined peaks are observed in 
rocking scans by rotating the sample-magnet assembly, 
which suggests that the magnetic flux lines are straight 
and parallel to the c-axis. 

The radial intensity distribution of the VL signal is dis- 
played in Fig. [2^. Two spectrometer configurations were 
used in order to maximize the signal-to-background ra- 
tio at different fields: collimation of 18 m with An = 10 A 
(configuration I) and collimation of 11m with An = SA 
(configuration II), where the sample-detector distance 
was equal to 20 m in both cases. The characteristic wave 
vector magnitude q of the VL is extracted by fitting the 
data to a Gaussian peak. This quantity is related to the 
magnetic field through the relation cr = (27r/q)^i?/$o> 
where $o is the flux quantum h/2e, and cr is a geomet- 
ric factor equal to \/3/2 or 1 for triangular or square 
VL, respectively. As can be seen from Fig. [^b, the mea- 
sured values of a are consistent with a triangular VL up 
to i? = 0.35 T (even though the determination of q using 
conflguration II has somewhat large errors because of the 
relatively poor momentum resolution). The geometry of 
the VL can also be determined from the azimuthal inten- 
sity distribution near the characteristic q. For B = 0.2 
and 0.25 T, the nearest diffraction peaks from the same 
type of domain are precisely 60° apart (Fig. , which is 
evidence that the triangular VL is undistorted. The fact 
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FIG. 1. Left: Neutron difTraction patterns at B = 0.2 and 
0.3 T and T = 2K, averaged over different sample rotation 
angles. Right: Rocking scans for two opposing Bragg reflec- 
tions (indicated by the arrows on the left), measured by rotat- 
ing the sample-magnet assembly with respect to the incident 
neutron beam. 
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FIG. 2. (a) Azimuthally-averaged signal amplitudes as a func- 
tion of momentum transfer q, measured at T = 2 K in differ- 
ent magnetic fields. Data obtained with two configurations 
(see text) are rescaled to allow a common vertical axis. The 
B — 0.35 T data are furthermore multiplied by a factor of 
5 for clarity. Vertical bars indicate the fitted value of the 
characteristic q. (b) The VL-structure-dependent coefficient 
a (see text) at different fields. Filled and open symbols corre- 
spond to configurations I and II, respectively, (c) Azimuthal- 
angle dependence of the signal measured in configuration I at 
T = 2 K for B = 0.2 and 0.25 T. Data are vertically offset for 
clarity. 
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FIG. 3. (a) Intensities measured at T = 2 K and B = 0.25 T 
in sector boxes surrounding the diffraction pealcs indicated by 
corresponding symbois in (b). Tlie data are offset for ciarity, 
and tlie same offset is used for refiections tliat are symmetric 
about tiie vertical plane, (b) Illustration of a tilt scan (see 
text) . Broken circles are tlie intersections of the Ewald sphere 
with the two-dimensional reciprocal plane at the end points 
of the scan. Potential diffraction spots are marked by crosses. 
The five diffraction spots marked by symbols have their Bragg 
conditions satisfied within the range of the scan, (c) VL in- 
tensity measured at T = 2 K as a function of magnetic field. 
Measurements in different configurations are normalized to 
the values at common magnetic fields. The dotted line is the 
expected field dependence in the London limit (Eq. the 
dashed line is the expected behavior after vortex-core correc- 
tion (see text), using Bc2 = 100 T (Ref. [H) and the phe- 
nomenological parameter s = 7; the solid line is a smoothing 
curve describing the data. The inset shows the same data 
plotted versus B~^''^, where the field dependence expected in 
the London limit is represented by the dashed line. 



that no intensity is observed between the peaks indicates 
that the VL domains are nearly fuhy aligned with the 
crystahographic axes. 

While it can be seen from Fig. [2^ that the signal 
decreases rapidly with increasing magnetic field, the 
azimuthally-averaged intensity does not provide the most 
accurate determination of the signal amplitude, because 
the Bragg conditions for different diffraction spots are not 
simultaneously satisfied. In order to determine the field 
dependence of the signal, we employed a procedure that 
is illustrated in Fig. ^jp: for each of the several magnetic 



fields measured in a given configuration, a scan is per- 
formed on the tilt angle of the sample-magnet assembly, 
which can be viewed as rotating the Ewald sphere about 
the r point with respect to the two-dimensional recipro- 
cal plane of the VL. By definition, both the Ewald sphere 
and the two-dimensional reciprocal plane contain the F 
point, and their intersection is a circle which expands as 
the tilt angle increases (arrows in Fig.|3b)- When this cir- 
cle goes through a diffraction spot, the intensity of that 
spot is maximized, and the different spots have their in- 
tensities maximized at different tilt angles, as can be seen 
from Fig.[2h- To measure the signal amplitude of the VL, 
we take the average of the maximal intensities of the five 
spots on the half of the reciprocal plain that corresponds 
to the tilting direction (for the B = 0.25 T example, the 
averaged value of the five maxima in Fig. [3^). In addition 
to configurations I and II, a third configuration III (col- 
limation of 18 m with An = 8 A and sample-detector dis- 
tance — 18 m) was used as a consistency check, and very 
good agreement is found among the data sets obtained 
in different configurations. The result is summarized in 
Fig.Et. 

The intensity Ihk of VL refiection (h, k) is given by2^ 
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where (p is the incident neutron fiux, /i is the neutron 
magnetic moment, V is the sample volume, and qhk is 
the magnitude of the wave vector of the refiection. F^k 
is the form factor, which is a Fourier component of the 
spatial variation of the magnetic field. In the London 
limit (i.e., without considering the effect of vortex cores), 
Fhk is related to the penetration depth A by 
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For B ^ Bel (the lower critical field), the second 
term in the denominator is dominant and this gives 
Ihk \/qhk B"^/^ (since for a given VL symme- 
try, qhk oc \/B). In type-II superconductors, the pres- 
ence of the normal-state vortex cores further reduces the 
VL signal at high fields, an effect that can be approx- 
imated by multiplying Eq.([2]) with the correction factor 
g-s S/Bc2^ where s is a phenomenological parameter on 
the order of unity and i?c2 is the upper critical field. 
As can be seen from Fig. [3]:, the measured intensity de- 
creases with increasing field more rapidly than cx _B^^/^ 
above B « 0.15 T (inset) and more rapidly than the esti- 
mated behavior after vortex-core correction above 0.25 T. 
The signal becomes undiscernibly small above 0.35 T, and 
measurements at higher fields (up to 4.0 T) do not indi- 
cate a reappearance of the signal (not shown) . 

The disappearance of the VL signal near B = 0.4 T can 
be explained by a mechanism^ first proposed for Bi2212: 
in an anisotropic superconductor, magnetic vortices are 
pinned over a characteristic length along the c-axis that 
is related to both the field strength and the anisotropy 
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of the penetration depth, 7 = (Ac/Aab)- At sufficiently 
high fields, this length decreases below the distance be- 
tween the adjacent CUO2 layers, and the VL decomposes 
into "vortex pancakes" that are two-dimensionally or- 
dered in each Cu02 layers but are uncorrelated along 
the c-axis. This substantially reduces the spatial varia- 
tion of the magnetic field and suppresses the diffraction 
signal. Such a dimensional crossover is expected to occur 
at B2D ~ *o/(7c)^ ~ 1.4 T for Hgf201, where c = 9.6 A 
is the c-axis lattice constant and 7 is '--^ 40 (Ref. [27I ). 
Both the estimated value of B2D and the measured value 
of 0.4 T are substantially lower than Bc2 — 100T,2i and 
than the magnetic field (B — 10.8 T) up to which the 
VL has been observed in the more isotropic compound 
YBCO," where i32D is expected to be much higher than in 
Hgl201. The difference between the estimated and mea- 
sured values of B2D might arise because 7 only serves as 
a phenomenological parameter in determining i?2D, and 
because there is considerable uncertainty in the measured 
value of 7 associated with doping^l and the experimen- 
tal method.— i^l"— We note that the width of the rock- 
ing curves in our measurement (Fig. [T|) suggests that the 
magnetic flux lines are straight over long distances. This 
finding is similar to the case of Bi2212 (Ref. 0) and can 
be reconciled with the dimensional-crossover mechanism 
for the VL decomposition if the flux lines are distorted 
only over short distances. In other words, as the flux 
lines break down into vortex pancakes, the magnetic fleld 
profile projected onto the basel plane dramatically loses 
its contrast, but the correlation between the averaged 
location of the flux lines (or stacks of vortex pancakes) 
remains intact. Finally, for the vortex pancakes that be- 
long to a single flux line to be misaligned along the c- 
axis, the pancakes in different Cu02 layers do not need 
to be individually pinned; instead, neighboring pancakes 
in the same CUO2 layer may be effectively pinned to- 
gether if the two-dimensional correlation length of the 
VL is larger than the size of each pancake and than the 
characteristic distance between neighboring pinning cen- 
ters in the same layer. Given the well-deflned diffraction 
patterns in Figs. [1] and [21;, it is likely that the vortex 
pancakes are pinned by a low density of pinning centers 
in this way. 

The temperature dependence of the VL signal was 
measured in two magnetic flelds B = 0.2 and 0.35 T. For 
each field, the measurement was performed at a fixed 
sample angle, and the intensities of multiple diffraction 
peaks were summed together in order to improve the 
counting statistics. Apart from a rescaling factor for the 
overall intensity, the temperature dependence measured 
in the two fields is identical (Fig. S]). Since B = 0.35 T 
is very close to the field at which the VL signal van- 
ishes, this suggests that the decrease of signal as Tc is 
approached (due to the increasing penetration depth, 
/ a A""*) and the disappearance of the signal at T = 2 K 
above B = 0.4 T have different causes. The similar tem- 
perature dependence in the two fields further supports 
the dimensional-crossover scenario of the VL decompo- 




Temperature (K) 



FIG. 4. Temperature dependence of diffraction intensity mea- 
sured at two different magnetic fields compared with /iSR re- 
sults for optimally-doped Hgl201,^ The value at T = 100 K 
is set to zero, and the data are normalized to the values at 
low temperature. The solid lines are empirical fits as guides 
to the eye (see text). 



sition, since B2D is related to the anisotropy 7 which is 
nearly independent of temperature.^^ 

The data in Fig. 2] can be described by the empiri- 
cal function / oc [1 - {T/T^f- The fitted value for 
Tc is 92 ± 5K, consistent with the Tc value of our sam- 
ple determined by magnetic susceptibility measurements 
(94 ± 1 K). The fitted value for a = 2.4 ± 0.3 is rather 
different from that of the "two-fluid" empirical expres- 
sion (a = 4) for conventional superconductors, which 
seems to describe early /j,SR results for YBCO.'^^ (In 
/iSR experiments, the measured muon relaxation rate 
is proportional to A^^, and therefore tr^ is expected 
to be proportional to the signal amplitude measured by 
SANS.) The data are also somewhat different from the 
more recent ^SR result for IIgl201)^ where cr^ can be 
described by the same function with a = 1.6 (Fig. [1]). 
These differences need to be further understood. We 
note that a discrepancy between SANS and /LtSR results 
regarding the temperature dependence has been previ- 
ously documented)^ and the temperature dependence in 
our data is consistent with other SANS measurements 
in the cuprates^ii^ For d-wave superconductors like the 
cuprates, a more quantitative description of the data near 
T = K may require a correction to account for nonlocal 
effectsi^ 



IV. DISCUSSION AND CONCLUSION 

Now we discuss the implication of our results for the 
VL physics in the cuprates. So far, well-deflned VLs 
at i? > IT have been observed only in YBCO and 
LSCO,^"^° and the results are drastically different: in 
YBCO,^ a hexagonal VL pattern is observed at low flelds, 
and a square pattern is observed only at B > 6.7 T, with 
the nearest-neighbor direction along the Cu-Cu diago- 
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nal of the unit cell. In contrast, a square pattern is 
already well developed at B ^ 0.4 T in LSCO,!'' and 
the nearest-neighbor direction is along the Cu-0 direc- 
tion, 45° from that in YBCO. A square VL has also been 
observed in the electron-doped cuprate Nd2-a;Ce^Cu04 
(NCCO, Ref. m down to B = 0.05 T. Our resuh for 
Hgl201 indicates that a unified picture may still exist for 
cuprates with high values of Tc (> 80 K), in which a tri- 
angular VL is favored at low fields {B < 0.5 T). We note 
that both Hgl201 and YBCO have non-body-centered 
crystal structures. Given that the formation of a square 
VL in YBCO is likely driven by the anisotropy of the su- 
perconducting order parameter;^ the square VLs at low 
fields in LSCO and NCCO do not seem to have the same 
origin, but may instead be related to their relatively low 
Tc values and/or to their common body-centered crys- 
tal structure. To further test the validity of this picture, 
measurements of the VL in Hgl201 and Bi2212 at higher 



fields (such as using local probes) are desirable. 

In conclusion, we have used SANS to observe the VL 
in Hgl201 near optimal doping. Our data suggest that 
a triangular VL at low fields {B < 0.5 T) is generic 
to cuprates with high Tc and non-body-centered struc- 
ture. The difficulty in using SANS to measure VLs at 
high fields in this anisotropic superconductor is likely 
due to a decomposition of magnetic flux lines into two- 
dimensional vortex pancakes. 
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